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Abstract 



A complete calculation of the finite temperature effective potential for 
the abelian Higgs model to the order e 4 , A 2 is presented and the result is 
expressed in terms of physical parameters defined at zero temperature. 
The absence of a linear term is verified explicitly to the given order and 
proven to survive to all orders. The first order phase transition has 
weakened in comparison with lower order calculation, which shows up in 
a considerable decrease of the surface tension. 



1 Introduction 

The electroweak phase transition fl], §, [3| has recently attracted much interest due 
to the hope to explain the observed baryon asymmetry of the universe within the 
minimal standard model ||, Much work has also been devoted to the phase 
transition in the abelian Higgs model, the simplest gauge theory with spontaneous 
symmetry breaking ||, [7], because it is believed to exhibit the main features of the 
electroweak phase transition. An important prerequisite to the understanding of the 
mechanism of the phase transition is the knowledge of the effective potential at finite 
temperature || §. It has already been calculated for the abelian Higgs model in 
|@, [1(J to the order e 3 , A 3 / 2 and in ]TI| to the order e 4 , A. In the last publication the 



assumption A ~ e 3 is made and the scalar masses m v and m x are counted as order 
A 1 / 2 , which can only be justified close to the critical temperature. Other approaches 



to the effective potential are found in [jE], jig] . 

In this paper, assuming formally A ~ e 2 , a complete calculation of the effective 
potential in the abelian Higgs model at finite temperature to the order e 4 , A 2 is 
presented. Coupling constants appearing in the denominator through infrared di- 



vergences |14| , p~5f| are taken into consideration and the full dependence on the Higgs 



field (f, the zero temperature vacuum expectation value v and T is kept. Therefore 



this calculation supplies the potential given in [TTJ] with A-corrections, which change 
the potential at T c and the surface tension significantly, in spite of their numerical 
smallness. 

The approach is based on the Dyson- Schwinger equation for the derivative dV/dip, 
i.e. the tadpole diagrams are summed H, [16| . In section 1 the principal method of 



the calculation, which makes also use of the gap equations for effective masses 0, 
is explained. 

The explicit formula for V((p, T) is given in section 2. It reproduces the results 
of II 10 and [11| when taking the appropriate limits. 



Section 3 addresses the problem of the linear term || [12|, [L7|, [18], [HJ. The 
expression for V to the order e 4 , A 2 shows that a nontrivial cancellation leads to 
lim^o dV/dip = 0, i.e. the absence of a linear term. It is proven that this feature 
survives to all orders of perturbation theory. 

To get rid of the arbitrary scale p,, which is introduced when regularizing the 
theory in the dimensional scheme, a transition to physical parameters is performed 
in section 4. However this finite renormalization at zero temperature proves to be 
not very important numerically. 

The numerical analysis in the last section concentrates on the surface tension 
PH , because this easily accessible parameter gives a first characteristic of the 
strength of the first order phase transition, which is important for the generation 
of the baryon asymmetry in the standard model. The surface tension is found to 
be generally smaller than the lower order calculation || suggests. This decrease 
is dramatic for large A, but in contrast to the e 4 , A-potential fill , the potential 
discussed here does not show a complete change to a second order phase transition 
in the considered domain of A. 
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2 Calculation of the Effective Potential using 
Dyson-Schwinger Equations 

2.1 Explanation of the Method in A(^ 4 -Theory 

After describing the calculation in some detail in simple A<£> 4 -theory the extension 
to the abelian Higgs model is shown to be straightforward. 
The euclidean Lagrangian has the form 

C=-\{d<p? + ^-\<p\ (1) 

where v = Xv 2 is counted as order A. Using the familiar zero temperature technique 
of Dyson-Schwinger equations [[EJ it is easy to obtain the relation 



--^z(V-V tree )=A + B = + 



Here the 3- vertex in the first term arises from the shift ip — > ip + (p. The two different 
sorts of blobs symbolize the full propagator and the full 3- vertex. 
The next step is to investigate the first term 



A = - X ^k 2 + m*+Tl(k) ' (3) 



where the tree level mass square mf ree = X0 — v is assumed to be of order A. For 
a calculation to order A 2 it is sufficient to know H(k ^ 0, k) and 11(0, ym) to order 
A and A 3 / 2 respectively. Here the need for different treatment of the contribution 
with zero Matsubara frequency can be understood by performing the substitution 
k — > ym in the integral. ( See 0, ITJ|, [T5[ for the correct way of counting the order 



of infrared contributions here and below. ) Therefore in the Dyson-Schwinger equa- 
tions for H(k) the vertex correction can be neglected : 



■n(A;) = -(n a (A;)+n 6 (A;))= + . (4) 
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With the definitions 



m§ = m^ ee + n a (o)| A3 / 2 + n 6 (o)| A , n 02 = n 6 (o)| A , (5) 
n 03 = n b (o)| A 3 /2 -n fe (o)| A , ni(fc) = n(fc) - n(o) , 

where the powers of A symbolize the accuracy to which a certain term has to be 
calculated, the following relation holds: 

m 2 tree + U{k) = m 2 + n 03 + Ui(k) + 0(A 2 ). (6) 

Of course H02 vanishes in A<^ 4 -theory, but keeping this term formally makes the 
extension to the abelian gauge theory with power counting rule A ~ e 2 more sim- 
ple. After substituting (|B|) into equation (^) the integrand is expanded neglecting 
contributions of order higher than A 2 , which results in the formula 



J \k 2 + rri\ (k 2 + m 2 ) 2 



dk 7-— 7; + 



3 {k 2 

n 02 n 02 + 1103 + 11! (A;) 1 



k 2 + m\ {k 2 + m 2 ) 2 ik 2 + m 2 ) 2 ) 

Here m 2 = m\\\ is the leading order mass term including the temperature correc- 
tion. Observing that in term B of equation (§) the vertex need not be corrected to 
obtain the full A 2 -result it becomes obvious that the third term in the right hand 
side of (0) together with B is equal to the derivative of the two-loop diagram 



£«>-£{ }■ 



with leading order mass corrections in the propagators. So the final expression for 
the potential is 

V = + f V V^ft + + V,. (9) 

A similar way of combining the different contributions to V has been considered in 
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2.2 Extension to the Abelian Higgs Model 



The abelian Higgs model in Landau gauge includes, from the topological point of 
view, exactly the same graphs, because it does also contain 3-and 4-vertices. Con- 
sider the euclidian Lagrangian 

C = --F^F^ - \D^\ 2 + z/|$| 2 - A|$| 4 , $ = ±={<p + ip + i X ). (10) 

Applying the formal power counting rule A ~ e 2 the main difference from the purely 
scalar model lies in the low order of the vector-scalar-scalar vertex. This manifests 
itself through the appearance of nonzero terms Ho2,t, ^-02, l with T, L referring to 
contributions to the transverse and longitudinal parts of the vector propagator || Q 
from the diagram in fig. 1. 

Fig. 1 

But this feature has already been accounted for in the last subsection by keeping 
n 2- Therefore the contributions to the potential can be classified analogously to the 
purely scalar case described above. This can also be shown starting directly from 
the Dyson-Schwinger equations of the abelian Higgs model. Of course at each step 
it is necessary to convince oneself that no contributions of order A 2 , e 4 are lost. The 
final result has exactly the same structure as the potential of the last subsection, 
but it involves more terms due to the particle content of the theory. 

V = V tree +V 1 + --- + Vs + V 3>11 + V 4 , u (11) 

Here 

V x = f dcp'Xcp'I(m x>3 ) , V 2 = f d<p'3\<p'I(m^) , (12) 

V 3 = f dip' (2 - 2e)eV J(m T>3 ) , V 4 = f dip'e 2 ip> 7(m L , 3 ) 
with the standard temperature integral 0] 

'M-fewh? (13) 

to be evaluated in n — 1 = 3 — 2e dimensions. The masses are defined in analogy 
with equation (||) by the appropriate zero momentum parts of the Dyson-Schwinger 
equations, which have to be iterated once to obtain the result to order e 3 , A 3 / 2 0. 
Corrections from the momentum dependent part of the Dyson-Schwinger equations 
are needed for the gauge boson only: 
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Finally V5 through V§ represent the diagrams of figure 2. 



Fig. 2 



Furthermore the contributions generated by the counterterm Lagrangian have to 
be added. 



3 Explicit Result to Order e 4 , A 2 

The calculations needed for the temperature dependent masses to order e 3 , A 3//2 can 
essentially be taken from ||. Differences arise because the loop with two propagators 
does only contribute in lowest order due to definition (|]). Also it is necessary to 
keep the e-dependence in leading order. The integrals in V5 through V§ have already 
been done in [11], [^J and wont be given here explicitly. After adding all the terms 
up, which proves to be rather laborious, the final result, using the MS-scheme, reads 



V{<p) 



„4 



(3 2 u A 



+ ^ + ? + ^K ln ^ 2 + l c + 3c2 
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+ 



16vr 2 
Ae 2 



6 8 64vr 2 



( In 2 r\ +w ^ + I In m+m *+ m x + 3 In 



(15) 



/3(2m+ mtp ) \ , fi 2 v\ 

3 ^ + 8tt 2 



(ln/i 2 /3 



2 c+f 



+ 



167T 2 

A 2 



(±ln/i 2 /? 2 + ±c-c 2 + ±-41n^ 



f3(m+mip+m x ) 3 ^m+m^) 

3 (m+m v ) 



+ 



(-§ ln/z 2 /3 2 + §c+c 2 + 2 + ln 



/?(m y +2m x ) 



16vr 2 

r A _J_ 

4 + 64^ 2 

4 



-3 In /3m,o + 2 In 



9 + 7^(10A 2 + 3e 4 )(c-|-ln/i 2 /3 2 ) + 



32vr 2 



rrr, 



'— In m, i 5 ( 2,n "^ rn "' 5 ) 
64tt 2 /3V (™+»v) 2 16tt 2 /5 



/Lf2 2 

IVI (T__ (2\ n P(-m+m<p+m x ) . j n 2m+m y \ 
T 2/92 \ 3 m+m^ / 



32vr 2 /? 2 
1 



(e 2 + X)m ip m x + e 2 m L (m„ + m x ) + 2\m(m x — mj) + e 2 m(m x + 2m tp ) 



12np 



m 3 + m 3 + 2m 3 + m|). 
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Here the lowest order masses are 



m 2 =3A<^ 2 + M 2 , m 2 =A<^ 2 + M 2 , M* = -v+ 4X + *f , (16) 



/l 2, 2 _2 2,„2 e 



12/3 2 

2 



Z Z Z Z, Zf Zf I 

m = e (p , m L = e + 



3/? 2 

and the constants c and c 2 arise from the temperature integral I(m) (see H) and 



the scalar two loop integral in (H) calculated in [25 



c= f + 21n47r-2 7 ^ 5.4076 , c 2 « 3.3025 . (17) 

Dropping terms of order e 4 (with power counting rule A ~ e 2 ) the effective potential 
from H is recovered. When changing the power counting to A ~ e 3 and m 2 , m 2 ~ A 
the result of Arnold and Espinosa |Tl| is obtained after dropping terms of order 



higher than e 4 . It should be noted that the term in (|IJ|) containing an explicit 
l/(^ 2 -factor does not show singular behaviour near ip = because of the logarithm, 
which decreases fast enough. 



4 Absence of a linear Term 

At first sight the contribution to the potential ( [L5[) proportional to m(m x + 2-171^) 
seems to produce a linear behaviour for small (p. But in this limit the logarithmic 
terms have to be expanded at the point (f = 0, which results in linear terms exactly 
cancelling the one mentioned above. In fact this feature, already discussed by several 
authors |9], 12], [17|, |18| starting from lower order calculations, can be shown to survive 



to all orders of small couplings perturbation theory : 

dV 

lim — — = to all orders in e and A. 

v^o dip 



This can be proven in the following way (compare the argument in [17|) : Global 
U(l)-symmetry implies 

1 8V n , n , . 

-^='4tf = 0) ■ OS) 

Obviously it suffices to show, that the self energy IT x (g 2 = 0) is finite for ip — > 0. 
Above the barrier temperature singularities can only arise from the transverse gauge 
boson propagator due to the temperature masses of x an d <p- Therefore diagrams of 
the kind shown in fig. 3 have to be investigated. Here the wavy lines symbolize full 
vector propagators and the blobs are vertices without internal vector lines, meaning 
the sum of all diagrams built from scalar propagators with the correct number of 
external vector lines. 
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Fig. 3 



If (p = 0, a gauge covariance argument, completely analogous to the zero temperature 
case, shows for the vertices with external vector lines only, that 

r 5...^ ~ N • • • \ k 2n\ for small |^| and fc? = 0, af3 . . . \w G {1, 2, 3}. (19) 

If <p 7^ 0, diagrams not covered by the previous gauge covariance argument because 
of explicit ^-factors at the vertices have to be added to T 2n . They however vanish 
not slower than <p 2 in the limit ip — > 0, which is clear from the fact that the unbroken 
theory has no vertices with an odd number of scalar lines. Therefore in the case 
| A; | ~ ip (p being the natural infrared cutoff introduced by the transverse vector mass 
ep>) the sum of both contributions, i.e. the complete vertex T 2n , can be counted as 
ip 2 when searching for small-<£> singularities. 

Consider the most dangerous lowest power of (p stemming from the maximal 
infrared divergence, which is obtained by setting k° = for all transverse vector 
propagators. It can be calculated by scaling the loop momenta according to k — > y<p. 
Counting the vector vertices as if 2 the following formula for the minimal overall 
power of p> is obtained (compare the argumentation in appendix A of fll5l ) : 

= 3L - 21 + 2(V - 2) + 2 . (20) 

Here L, I and V denote the number of vector loops, vector propagators and full 
vertices, symbolized by blobs in fig. 3, respectively. The last term +2 follows from a 
closer look at the contribution of the full vertices with external x-lines: Fig. 3 shows 
examples of the two different structures to be investigated. If there are two such 
" X- vertices" , each will contribute a factor ip, resulting in the correction +2. If there 
is only one with two external x-hnes, it may have no explicit ^-factor. In the latter 
case however the contribution of the vector vertices 2(V — 2) has to be replaced by 
2{V — 1), which again corresponds to a correction +2. Therefore equation ( p0|) is 
valid in the general case. Now using the well-known formula V + L — I = 1 it follows 
immediately that 

?V = L>0, (21) 

or equivalently: There is no divergence for (p — > 0. If not all of the vector propagators 
are infrared divergent, the vertices connected by "heavy" lines may be formally 
fused. Now repetition of the above argument leads again to the desired result thus 
completing the proof. 
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5 Transition to physical Parameters 

To get rid of the arbitrary scale Jx the potential is rewritten in terms of physical 
parameters defined at zero temperature. Such parameters are the Higgs and vector 
masses and the vacuum expectation value of the Higgs field. To stay closer to the 
previous notation they can be expressed through new coupling constants A and e, 
defined by 

dV 



2 _o\2 2 _-2 2 

m V ,phys — LAV i m phys ~ e V ) ~ 



= 0. (22) 

^Pphys ^ 



Effectively a finite renormalization of the form 

<P 2 = ¥> 2 ph y S (l + c), A = A + 5A, (23) 

'^phys 



e 2 = e 2 + <5e 2 u = Xv 2 + 5u v 



has to be performed. The physical (^-propagator is 

1 



(24) 



(l + c)(g 2 -m 2 -lW))' 
therefore the on-shell definitions of the new parameters follow from 

c = Atonal (25) 

'-'y * ip.pnys 

< + n,K, p% ,) = m^ ht „ m> + Re n(^„) = m \ h „ 



together with the last equality in (E3J). The zero temperature effective potential to 
the order A 2 , e 4 needed here is given by 



A a ™t 



V = -_^ + - V9 4 2L | + l n 4 (26) 

2^ ^ 4 y 64tt 2 V 2 m x/ v ; 

64tt 2 l2+ m <; 647r 2 U + ln m 2 J- 

Contributions to the self energy corrections IL^g 2 ) and IT(g 2 ) come from the usual 
zero temperature one loop diagrams. In view of the principal features of the potential 
considered in this paper the numerical effect of the performed finite renormalization 
is not very important (see fig. 5 in the last section). Therefore the complete formula 
for V, which is easy to obtain, is not given here explicitly. 
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6 Numerical Results and Discussion 



In the previous sections a complete calculation of the finite temperature effective 
potential to the order e 4 , A 2 has been performed, including the transition to physical 
parameters defined at zero temperature. The gauge coupling is chosen to be e = 0.3 
and the influence of A, which corresponds to the Higgs mass, is investigated. At 
first sight the calculated potential seems to ensure a first order phase transition in 
a wide range of the parameter A, but reliability of perturbation theory has to be 
questioned. This becomes obvious from fig. 4, where different approximations of 
the effective potential at their respective critical temperatures are shown. (Here 
and below the dimensionful quantities are given in units of v and its powers.) Some 
insight can be gained from a comparison with the results to order e 3 , A 3//2 obtained 
in ||, [l0| . A reasonable physical quantity to be calculated from both potentials is 
the surface tension [T2I], ^TJ , which can be seen as a measure of the strength of the 
phase transition : 



a = j o dip^2V(<p,T c ), (27) 

where ip + is the position of the second degenerate minimum of V and the potential 
is normalized to ensure V(ip = 0) = 0. The numerical results are shown in fig. 5. It 
includes besides the surface tension from potentials to the order e 3 , A 3 / 2 [[| and to 
the order e 4 , A 2 also the results calculated from a potential to order e 4 , A [11|, where 



according to the power counting rule A ~ m 2 ~ m 2 ~ e 3 all terms of order higher 
than e 4 have been neglected. Obviously the shift introduced by the transition to 
zero temperature physical parameters is not important for the present discussion. 

The fact that perturbation theory is not reliable for large A, already stressed 
in ||, can be clearly read off from fig. 5 (a changes by an order of magnitude 
when adding the last term in the perturbation series). Somewhat surprisingly the 
perturbation series cannot be trusted too much for small A either. Here a decreases 
by a factor of ~ 2 at least. Two different higher order terms are mainly responsible 
this for this change. In the region of extremely small A it is essentially the e 4 y? 4 term 
(see equation flTSP), which cannot be viewed as a small correction to the tree level 
term \(p 4 /4. This enhancement of the y9 4 -term is a temperature effect and therefore 
cannot be removed by zero temperature renormalization. For m 2 /m 2 ~ 0.4 the 
logarithmic mass dependence of the e 4 ^ 2 term seems to be more important. The 
great influence of this numerically small contribution can be understood by recalling 
that at the critical temperature the leading order <£> 2 -terms essentially cancel and 
that a ^-dependence in a coefficient of <p 2 /f3 2 cannot be absorbed in a correction 
ofT c . 

It is interesting to compare the above discussion with another method to investi- 
gate the reliability of the perturbation series: Applying the ^-conditions introduced 
in [H] in the more restrictive form of |15|] (i.e. using the higher order expression 
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for the masses) the region of reliability m 2 /m ~ 0.1 is obtained for £=2. For 
the largest permissible Higgs mass fig. 5 suggests Aa/cr = 0.57 which signals the 
breakdown of the perturbation series. This error does not decrease for smaller Higgs 
masses, in contrast to the standard model calculation in lower order of JTH], due to 
the shift of A through an e 4 -term discussed above. However, this enhancement of 
the (/^-contribution, which is invisible in the ^-conditions, does not threaten the first 
order of the phase transition. Therefore perturbation theory seems to ensure a first 
order phase transition for small A at least, in spite of the still unknown exact value 
of the surface tension. 

The potential to the order e 4 , A from [11] does not give rise to a first order 



phase transition for A ~ 0.01. The value of a calculated using this potential differs 
significantly from the result presented here for A ~ 0.007. This is partially due to 
contributions of the form e A f 2 ln(m + m v ) and the like, already mentioned above. 
Counting as higher order correction results in the In m-contributions found in 



the potential from [|TT], which is obviously a significant change for small ip. 

Fig. 6 shows the vacuum expectation value in the asymmetric phase at T c . This 
parameter does not reflect the dramatic change of the surface tension by higher 
order corrections. 

It would be interesting to extend this approach to the standard model, which 
seems to be straightforward, and to try to estimate the influence of expected cor- 
rections beyond e 4 , A 2 . 

I am most grateful to W. Buchmiiller, who suggested this investigation, for 
continuous support and encouragement. Also many helpful discussions with 
D. Bodeker, Z. Fodor, T. Helbig and H. Kohrs have to be acknowledged. 



Fig. 4 Different approximations of the effective potential plotted at their respective 
critical temperatures with A=0.01 ( the e 4 , A-potential is a result of [|TTJ ) 

Fig. 5 Dependence of the surface tensions calculated from the different potentials on 
the zero temperature mass relation m 2 /m 2 = 2A/e 2 with e = 0.3 

Fig. 6 Dependence of the position of the second minimum of different potentials on 
the zero temperature mass relation m 2 /m 2 = 2A/e 2 with e = 0.3 
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